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We present a general method for analysing novel computational 
substrates to determine which of their parameters can be ma- 
nipulated to exhibit the complete set of 2-input boolean logical 
operations. We demonstrate this approach with an NMR-based 
case study, showing which NMR parameters can be used to per- 
form boolean logic. 
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1 INTRODUCTION 

There is much work on in materia computing: exploiting unconventional 
material substrates and their dynamical properties to perform computation, 
and examining their computational properties and abilities. One of the more 
advanced is optical computing ll33l[34l . Other, more exotic, substrates include 
nuclear spins in NMR experiments |[l6l[l7l|24l|6l, Hquid crystals llT2l[T3lfT4ll . 
conductive media Ii21u22j . reaction-diffusion chemical systems 1791 l23l l27l 
[321 . DNA EH), and even slime moulds Il2ll3l1. 

The aim is to allow the material to do what comes naturally, under control 
of physical laws, and to exploit this natural dynamics as an efficient form of 
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class members 

T, F 

1 A, NOT A, B, NOT B 

2 A AND B; A NAND B; NOT A AND B, A AND NOT B,... 

3 A XOR B, NOT {A XOR B) 



TABLE 1 

The four equivalence classes of the 16 boolean 2-input logic gates. These correspond 
to equivalence classes under permutations of inputs and negation of the inputs and/or 
outputs f9i, and to the different kinds of canalising functions 1.101 . 



computation l|29ll30l[T8l . The question arises: given some novel substrate, 
how to analyse its properties to determine whether it is a suitable medium for 
computation, and, if so, how it can be manipulated to perform computation in 
a manner best suited to that substrate. 

Here we tackle a part of that problem: analysing a novel substrate to 
determine how it can be used to perform boolean logical operations. We 
present the design principles through a case study of using bulk nuclear spins, 
in the context of NMR (nuclear magnetic resonance) experiments. 

2 THEORY AND BACKGROUND 
2.1 Boolean functions and logic gates 

The NAND and NOR logic gates are both universal, in that either can be 
used to construct all the other boolean 2-input logic gates. However, to do 
so, several such gates may be required. In an implementation, it is often 
more important to minimise the circuitry (the number of actual gates) than to 
minimise the number of types of gates. 

Logic gates can be classified in terms of their symmetry properties ll28l[m . 
In terms of minimising circuitry, gates can be classified in terms of equivalence 
classes on permutations (rewiring) of inputs, and on negation of inputs and/or 
outputs (i.e. adding inverters, important in cases where they are significantly 
cheaper to implement than binary gates |^). There are four such equivalence 
classes for the 16 boolean 2-input logic gates (Tables [T] and |2]l. 

The same classification is found by considering canalising functions ifTOl . 
These are functions where the output is independent of one (or more) of 
the inputs, for one (or more) input values. Class consists of the constant 
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TABLE 2 

Truth tables of examples of the four classes of boolean logic gates, for each possible 
input value of A and B. Class (example T) is independent of either input value. 
Class 1 (example B) is independent of one of the input values (here, A) for all possible 
inputs. Class 2 (example A NAND B) is independent of one of the input values when 
the other has a particular value (here, it it independent of the value of B when A has 
the value 0). Class 3 (example A XOR B) depends on both input values for all inputs. 



Class 




, B 1 NAND 4 XOR 



□ 



FIGURE 1 

Logic-gates patterns and their different symmetries as visible from the two-dimensional 
representations of the canalising functions gate classes 0, 1, 2 and 3. 



functions: the output is completely independent of the inputs. Class 1 consists 
of the strongly canalising functions: the output is independent of one of the 
inputs (for example, the function B is independent of the value of input A). 
Class 2 consists of the weakly canalising functions: the output is independent 
of one of the inputs when the other input has a specific value (for example, the 
function A NAND B is independent of the value of input B when the value of 
input A is false). Class 3 consists of the non-canalising functions: the output 
is determined by both inputs. 

We use these symmetries/canalising properties of the gates (Figure to 
look for analogous properties exhibited by the substrate that indicates a natural 
implementation route. 
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2.2 Nuclear spin dynamics 

The most common application of NMR spectroscopy is that as a tool for 
structure elucidation of condensed matter in general, and of molecules in 
the liquid state in particular. The quantised nature of nuclear spin can also 
be assigned to the notion of a qubit in quantum computations 1261 . Nuclear 
spins have also been used as a substrate to implement classical computation 
paradigms such as binary or continuous logic gates and circuits |24, 6|. 

A major reason for these computational applications to choose nuclear 
magnetic spins as the implementation platform is the rigour by which nuclear 
spin dynamics are described by quantum mechanics. In addition, the accuracy 
by which the macroscopically measured signal in an NMR experiment can be 
related to the underlying microscopic spin dynamics, using a density matrix 
approach, is nearly unrivalled by any other spectroscopic technique. In order 
to develop the formal relationship between the properties of logic gates and 
the quantum mechanical expressions describing the spin dynamics in a NMR 
experiment, we give a short summary of the necessary formalisms. For more 
details see textbooks on NMR |20|. Here we consider only uncoupled spins 
5* = 1/2 in the liquid state. 

The time evolution of an initial spin state vector (ta)) is described as 

\^P{t)) - Uit,ta)\^ita)) (1) 

where U {t, ta) is the time propagator describing the spin dynamics at every 
given point in time. The general orientation of the spin vector \ tfj) is ll20l 

where 6^ and 0s are the polar and azimuth angles. Its three-dimensional 
representation is given by 

{■^\Ix\ip) \ , / sin COS \ 
(V'l/ylV') - 2 Sine, sin (3) 
{i^\h\i>) J \ cos9s J 

describing a general Cartesian vector orientation (see Figure |2ji, where I^, ly 
and Iz are basis spin operators ||251 . 

The time propagation of the macroscopic NMR signal is described by the 
density matrix p{t) = ji/i {t)) {tp {t) \ as 

p{t) = U{t,ta)p{ta)UHt,ta) (4) 
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FIGURE 2 

Spin state vector orientation of j ^) in its Cartesian representation 



where the bar signals the ensemble average. The thermal equilibrium density 
matrix is ll20l 

p^.'i- ^I + ^Ab/. (5) 

where Ab = H'^Bq/ [k-QT) defines a Boltzmann factor, scaling the separation 
of spin energy levels, that can be interpreted as a spin polarisation along the 
z-axis. 

The NMR system Hamiltonian and generator of the time evolution propa- 
gators for isolated, uncoupled spins is 

n^ncs + n,{ (6) 

where Ties and Hrf are the chemical shielding and radio frequency (r.f.) 
Hamiltonian, respectively. The explicit representation of the rf. Hamiltonian 
is particularly simple in the rotating reference frame (RRF) as 

Hrf = ujph + Hp (4 COS (f>p + ly sin (j>p) (7) 

where Wp is a frequency offset (relative to the frequency of RRF), and Kp 
and (j)p are the amplitude and phase of the pulse. Figure[3ja) summaiises the 
behaviour of the magnetisation vector M (Eqs. ([T]i-(|23]l) under the influence 
of a r.f. pulse; Figure[3jb) depicts the relevant parameters. Pulses applied to 
spins resonating at the rotation frequency of the RRF are described by the 
system Hamiltonian 

n = "Hrf = Kp (4 cos 0p + ly sin 0p) (8) 
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r.f. pulse delay signal 




FIGURE 3 

NMR single-pulse experiment: (a) rotation of the magnetisation vector M by a r.f. 
pulse i?0p with(^p — ^ and/3 — KpTp = |-; (b) the relevant parameters describing 
the single-pulse NMR experiment (pulse amplitude Kp, pulse duration Tp, pulse phase 
(l)p, pulse frequency ujp, free evolution delay Td and acquisition phase <^a). 



and the overall system Hamiltonian for uncoupled isolated spins S = 1/2 is 

= 'Hcs + 'Hrf : during pulse 
^ = ^cs : during free evolution (9) 

including a chemical shielding offset Hcs from the RRF frequency. 
The general form of the time evolution propagator in NMR is 

U{h,ta) = Tcxpj-i^ ''H(t)dt| (10) 

where T is the Dyson time-ordering operator ||251 . 

Dependent on the symmetry of the Hamiltonian, different (simpler) expres- 
sions for the propagator can be formulated 

H{t)=H ^ U{h,ta)=exp{^m{h~ta)}(n) 
time independent 
[Hit'),H{t")]^0 U{tb,ta)^exp{-iJl\\t'H{t')}(l2) 

self — commuting 
{ta +t) ^ -n (tf, ~-t) U (tb, to) = 1 ; anti-symmetric (13) 

W (4, ta) =±1; cyclic (14) 

In the case of a hard pulse, a so-called S pulse is usually a valid description 
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of the time propagator 



U {tt, to) = CXp {-iHrfTp} ; Tp = tb-ta 

= exp{-iKprp (/j; cos(?!)p + sin(/)p)} 

= R^JP);P^KpTp (15) 

which is equivalent to the rotation operation R^^ (/3). It generates a rotation 
by an angle /3 about an axis in the zy-plane with azimuth angle (pp. For a 
sequence of 6 pulses the propagator can be factorised as 

= (1^2) Rcj,^! ; f3l = KplTpi; /?2 = Kp2Tp2 (16) 

The explicit form of these propagators is equivalent to rotation operators in 
two dimensions and is 1201 



-i/3(/^ cos0+/ysin0) _ f COS ^ — isln^e 
V -isinfe+''^ cosf 



-R0(/3) = e- 

= R,{.p)RAl3)RA-^) (20) 



After a pulse R^^ with phase (j)p ~ ^ and flip angle /3 = f the thermal 
equilibrium density matrix pf^ is transformed according to 

p. (21) 

yielding spin polarisation along the +x direction. 

A r.f. pulse R^^ applied to an initial state p*^*) generates a spin density 
matrix that can be expressed in terms of its final polar coordinates p (0s, Og) = 
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(V'l (Eq. ([2]l) or as a function of the pulse parameters rotating an initial 
state {(ftp, (3). For example when starting from the final state is 



a\ _ ^Ti , -^^ f -isin/3c "^p 



and the relation p {(j)s, Os) ~ p*-^-* {4>p, (3) holds for (0s, = ((/)p — §,/?)• 

Equivalently to the calculation of the three-dimensional spin-vector ori- 
entation in Eq. ^ from the two-dimensional spin state lip) (Eq. (|2|), the 
three-dimensional magnetisation vector M can then be calculated from any 
given two-dimensional density matrix p as 





, / Tr{p4} 




= Tr{p/,} 


ih) ) 


' V Tr{p4} 



M = (4) = Tx{plA (23) 



In the case of single isolated spins this also implies the equivalence in interpre- 
tation of the spin vector orientation in Figure [2] and the magnetisation vector 
M. 

In general a NMR pulse sequence can consist of any number and parametri- 
sation of r.f. pulses and free evolution delays. This provides for numerous 
options for the implementation of e.g. logic gates. The most useful parametri- 
sation in the following is the one describing the NMR signal as a function of 
the pulse sequence, p^*' (</)p, 



3 LOGIC GATES FROM NMR PARAMETERS 
3.1 Classification of NMR parameters 

The most common parameters describing a NMR experiment, open for use as 
control parameters to implement logic operations, are {^p, 0p, Tp, Wp, Td, ^aS'- 
pulse amplitude, pulse phase, pulse duration, pulse frequency, free evolution 
delay, and receiver phase, respectively (see Figure|3]l. 

A heuristic classification rule has been given elsewhere [24 1, separating 
logic gates based on the relation between the experimental NMR parameters 
chosen: 

1 . if the effect of the first parameter cannot be compensated for by set- 
ting the second parameter (canalising input value), then AND, >, < 
NOR, OR, <, >, NAND gates can be constructed (asymmetric gates); 
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A NAND B 
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TABLE 3 

Number of canalising input values for the four gate classes. Demonstrated for one 
representative member of every class: T, B, A NAND B, A XOR B. The respective 
truth tables for these four gates are shown in Table|2] 



2. if the effect of the first parameter can be compensated by a setting of 
the second parameter, then XOR and XNOR gates can be constructed 
(symmetric gates). 

This classification of logic gates seamlessly integrates in a generalised descrip- 
tion by using the concept of canalising functions lITOl and equivalence classes 
121 . Fixing one of the inputs of a logic gate {A or B) is a canalising input 
value if changing the other (unfixed) input does not alter the gate output. The 
number of all possible canalising input values for a given logic gate for all 
possible inputs states is used as a criterion to assign logic gates to one of four 
classes (see Table [3]l: 

class two canalising input values for either of the two inputs A and B 
class 1 two canalising input values for exactly one of the inputs Aor B 
class 2 one canalising input values for either of the two inputs A and B 
class 3 no canalising input values for either inputs A and B 

These four classes also show unique symmetry patterns in their two- 
dimensional representations (see Figure [TJ. Rows and columns of equal 
coloured squares signal the presence of a canalising input. These patterns can 
be used to map functions f{xl, x2) to corresponding canalising functions and, 
therefore, the logic gate(s) they can naturally implement are readily identified. 

According to our previous definitions [24 J a strong parameter would gener- 
ate a canalising function for at least one of its states (class 0, 1 and 2), because 
altering the other parameter causes no change in the NMR output. On the 
other hand, for a weak parameter the output would always change (class 3). 
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In the following we demonstrate how the abstract concept of canalising 
functions relates to the symmetries and commutation relations of the NMR 
Hamiltonian, the time evolution propagator, and the experimental output of a 
NMR experiment. 

3.2 Canalising input values in NMR context 

Here we consider a system Hamiltonian that consists only of the Hii term. 
Hence, we are in the single spin, strong pulse, on-resonant regime, and for 
simplicity we only consider the experimental parameters Kp,(j>p, Tp (r.f. pulse 
amplitude, phase and duration) for logic gate generation. 

The NMR r.f. pulse propagator R^^ = R, (0p) R^ R, (-0p) 
(Eq. ( pOj l) is the general rotation operator about an axis in the traverse xy- 
plane of the RRF. It therefore commutes with a particular spin state (/^p) in 
the transverse plane for a given value (pp (e.g. for (j>p = 0; [Ro{f3), Ix] = 
independent of the value of f3). However, i?^^ (/3) will never commute with 
the thermal equilibrium state pf^ (~ /^) (Eq. (|5]l), which is perpendicular 
to the aiy-plane. For example, a i?o (f ) pulse applied to z-magnetisation 
(^ Iz) will flip it to the —y direction, and therefore change the spin state 
and its orientation. A subsequent pulse apphed using (f ) will leave it 
unaltered, since the system is in the —ly spin state, an eigenstate of R~^ (f ) 
and, therefore, a canalising input is generated in this second step. One can 
postulate now that a canalising input is generated by an r.f. pulse if: 

1. {(3) leaves the system state unaltered (eigenstate; p] — 0), 

2. — 1, the unity operator. 

In the following this behaviour and the possibility of generating canalising 
input values is analysed. First we consider single-pulse experiments assuming 
different initial system-state preparations p^*-* . Second we consider experi- 
ments composed of more than one pulse. 

Single-pulse gates 

The signal as detected during a NMR experiment is not the effect of the full 
spin magnetisation vector M as given in Eq. ( |23| l but only its projection into 
the xy-plane [20 1. Hence, the physically meaningful quantities that can be 
measured by NMR are the magnetisation vector components = (Ix) 
and My ~ (ly), describing the orientational distribution of this projection, 

and the magnitude Al^y = \J {I-x)^ + {ly)^ describing the total amount of 
magnetisation present in the a;y-plane. The functional structure of M^, My 
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and Mxy is determined by the r.f. pulse R^^ (/?) and the initial spin state p^*^ 
to which the r.f. pulse is applied. Mx, My and M^y are therefore functions of 
pulse phase (f>p, pulse flip angle /3 — KpTp and the initial direction of the spin 
magnetisation. 

The spin magnetisation can assume every possible orientation, while the 
rf pulse and therefore the rotation axis of the spin magnetisation is restricted 
to the xy-plane. It is instructive to examine scenarios where the initial spin 
state p*^*^ is either perpendicular to the a;y-plane and never commutes with 
the r.f. pulse operator R^^ (/?), or where it is coplanar to it and therefore can 
commute. 

Starting from thermal equilibrium state The most simple and natural ini- 
tial spin state is the one the system assumes at thermal equilibrium p*^*^ = pl'^ 
(Eq. (|5]l). Here the spin polarisation is pointing along the z-axis, perpendicular 
to the xy-plane and R^_^ (/3) never commutes with it. p^'^ transforms under the 
influence of a pulse R^^ as 

p(^)(0p,/3) = R4>MpT,Rlm (24) 

After this single pulse the spin magnetisation vector M (Eq. ( |23] l) is a function 
of pulse flip angle /3 and pulse phase (pp 

\ / Tr{p(^)((/.p,/3)/,} \ / sin sin /3 \ 

My = Tr{pW(0p,/3)/,} =^ -cos^pSin/3 (25) 
y V Tr{pW(0p,/?)/4 / V cos/3 / 

and the NMR measurable quantities AI^ {4>P7 My {(f>p, (3) and M^y {(f>p, (i) — 
^ I sin j5\ are readily determined. These quantities are functions of the two 
pulse parameters 0p and /3 just like the binary logic gates are functions of 
the two inputs A and B (Table |2|l. An instructive way to analyse M^ (0p, 
My {(j)p, /3) and Af^j, {(j>p, /3), as to which logic gates can be implemented 
by them, is by examining their representation as two-dimensional contour 
plots. These are shown in Figure]?] A comparison with the two-dimensional 
sketches of 2-input logic gates in Figure]!] immediately reveals agreements 
and disagreements in symmetries. 

In principle the functions Mr^ {(j)p, f3). My (0p, /3) and M^j, {(j>p, f3) are 
continuous in (j)p and /3 whereas logic gates are boolean functions of A and 
B that can only assume the discrete values of {0, 1}. In order to map the 
continuous functions in Figure]4]to the four possible discrete sets (A,B) of a 
logic gate one has to find four discrete pairs {(pp, f3) at which to evaluate the 
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FIGURE 4 

Contour plots of the NMR detectable quantities (a) Af^ = (J^), (h)My = (ly) 

and (c) Afj^y — \f(Jx) + {ly) (Eq. \25\ ) as a function of the pulse parameters 
(f>p and /3. The initial spin polarisation has been along the z-axis ip"^) and Ab = 1 
has been assumed. The ranges for (j)^ and /3 have been chosen such that the only 
symmetry operation necessary to generate a full, infinite plot are horizontal and vertical 
translations. 



continuous functions and where the symmetry pattern of the desired logic gate 
results. 

For example from Figure |4]^ a) depicting M^, one can directly identify 
the symmetry pattern corresponding to a class 3 gate in Figure [T| In or- 
der to implement the class 3 (XOR) gate (Table |2| one needs to define four 
discrete value sets (0p = (/ip ,/3 = /3^) (with A,B e {0, 1}). Inspection of 
Figure |4]^a) immediately suggests the positive and negative extrema as possi- 
ble candidates: Selecting the pair (0p /3^=°) = (f,-f) the detected 
output is (0;^=" = f ,/3^=° = -f ) = -0.25. In the same fashion, 

(f , f) = 0.25, ,-f) = 0.25 and , f ) = -0.25 are 

calculated for the remaining XOR gate input configurations. The values of 

(0p, My ((/>p, j3) and M^y (0p, (3) in Figure|4]are calculated assuming 
Ab = 1. Mapping these values to the final boolean values {0, 1} is achieved 
by scaling them by a factor of 4. 

A class 3 gate is characterised by the absence of any canalising input 
(Table p|. A change in (j)^ or (3^ will therefore always change the value 

[<t>p, f3^)- In order to change one of the parameters (p^ or and achieve 
the same result in {'Ppi P^) ons can take advantage of the periodicity 
of the trigonometric functions in Eq. ( |25| ) and in this way ( f 7 — f ) = 
-f ) = -0.25 and (f , f ) = f ) = 0.25 can be im- 

plemented. This gate is evaluating (p^ at values separated by 2tt and the 
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TABLE 4 

Some selected values 0p and /3 to generate one representative gate of every gate class 
using Mx ((^p, /?) as shown in Figure[4|a). 



gate is therefore independent of a change in (p^ no matter which of its two 
permitted values assumes. A change in f3^, however, will always change 
Mx (0p , /3^). This behaviour corresponds to a class 1 gate with two canal- 
ising input values for one parameter. In order to implement class 2 gates a 
comparison of the logic gate patterns and the contour plot for ((/)p , 
suggests a restriction of the parameter range to smaller intervals. For example, 
one may choose the interval (pp g [t, ^] and f3 E [O, f ] to implement a 
NAND gate representing a class 2 (NAND) gate (see Table |4]i. The implemen- 
tation of a class gate requires two canalising input values for both 0p and 

. As in the class 1 (B) gate scenario the periodicity of {4>p7 can be 
used to achieve canalising inputs. However, this time both parameters have to 
undergo the 27r value changes. Alternatively one can use pulses corresponding 
to rotation operators i?^^ {/3) = ±1 for generating canalising inputs for any 
parameter configuration of <j)p or (3^ . 

Having demonstrated how to implement one member of every gate class 
it is now a trivial task to generate implementations of all the members of a 
given gate class. By using the symmetry operations of permuting the input 
assignment of A and B to (pp and /3, inverting the input or the output, the 
remaining gates are obtained directly. 

This shows that all 2-input logic gates can be implemented using the 
functions (0p, /3). AIj. (0p, (3) and My {4>p, /3) are both products of two 
linear trigonometric functions, differing only by a phase shift of | in the factor 
depending on (pp (sin (0p + f ) = cos^p). This phase shift can be seen in 
Figure|4|a) and (b) as a horizontal shift by ^, otherwise both plots are identical. 
Their overall symmetry pattern is that of a class 3 gate. Hence, everything said 
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about Mx {(j)p, /?) implementations equally holds for My {<j>p, j3). 

However the same can not be said about the magnetisation magnitude 
M.,y ((/.p,/?) (Figure Qc)). M,y (d^p,l3) = M,y (/3) = ^ |sin/3| is only a 
function of the pulse flip angle /3 and, therefore, displays a contour plot with a 
symmetry corresponding to a class 1 gate. 

The behaviour of the functions under permutation of the variables (p-p and (3 
and/or their inversion distinguishes ((/>p, My {(j)p, /?) and M^y {4>p, (3) 
from each other AI^ (0p, /?) is invariant to permutation of the variables while 
Mxy {(t>p, (3) is not. In contrast, M^ {(f>p, f3) changes under inversion of its 
variables (f)p — —(j)p or /3 — >■ —(3, while Al^y {ipp, (3) is invariant to inversion 
of (f>p — > — 0p. My {(j)p, (3) shows the same behaviour as M^y {<j>p, (3) under 
these transformations. However, a simple shift in (/)p = ^ gives Mr {(j)p, (3) = 
My (</)p + ^ , /3) , an operation which is not possible for M^y {(f>p, (3) since it is 
independent of (pp. The differences in behaviour arise because for M^ (0p, f3) 
and My (0p, (3) both variables 0p and /3 are arguments of products of trigono- 
metric functions, which generates function values in the range of [—1,1], 
whereas M^y {4>p, (3) has trigonometric function values in [—1, 1] only for the 
variable /3. 

In principle, the task of implementing all 2-input logic gates by NMR 
spectroscopy is already accomplished by using only the simplest of all NMR 
experiments: starting from thermal equilibrium and only using a single r.f. 
pulse. This scenario can therefore serve as a Universal Logic Module (ULM) 
||35]| . There are, however, good reasons why one needs to explore other 
starting conditions, and more complicated NMR pulse sequences. NMR is 
uniquely suitable as a single testbed for the implementation of classical as well 
as quantum computations. In common formulations of quantum algorithms 
the initial state of the computation is not the thermal equilibrium state but a 
superposition state 1 16 |. For valid comparisons between classical and quantum 
algorithm NMR implementations the same initial (superposition) state should 
be used. Further, it may be desirable to construct more extended circuitry 
than just a single logic gate ITSl . Then gate implementations that can deal 
with initial states other than p'i^ are attractive as an efficient means of taking 
advantage of the output of a logic gate without having to restore the initial 
state before the next logic gate in a circuit can be executed. 

Starting from superposition state A r.f. pulse i?f (f ) applied to p''^ gen- 
erates the superposition state p^ = \Tl + which corresponds to spin 
polarisation pointing along the +a;-axis. Here we take p^. as the initial super- 



14 




0.5 1 
p[rad]/TT 




0.5 1 
</)p[rad]/7T 




FIGURE 5 

Contour plots of the NMR detectable quantities (a) Alx 



(h)My = (ly) 



2 2 

and (c) Mxy ~ \/ (Ix) + (ly) (Eq. \2S[ ) as a function of the pulse parameters 
0p and /3. The initial spin polarisation has been along the a;-axis (px) and Ab = 1 
has been assumed. The ranges for and /3 have been chosen such that the only 
symmetry operation necessary to generate a full, infinite plot are horizontal and vertical 
translations. 



position state. The commutator between i?^^ and pj, is 

[R^^{l3),p^] = -XbIz sin (j)p sin ^ (26) 

and is zero for values 0p — nir. These angles correspond to rotations around 
the positive and negative x-axis (the commutator is also zero for j3 — 2nn, but 
these angles represent trivial 27r rotations), transforms under the influence 
of a r.f. pulse R^^ as 

p(^)(0p,^) = R^^{P)p,Rl^{l3) (27) 

The r.f. pulse generates a magnetisation vector M (Eq. ( |23] l) as a function of 
pulse flip angle /3 and pulse phase 0p according to 

M, \ / Tr{p(-)(0p,/3)4} \ , / l-2sin2 0psin2§ \ 

= Tr{p(-)(0p,/3)/,} =^ sin20psin2f (28) 
/ V Tr{p(-)(0p,/3)/4 / V -sin0pSin/3 / 

All the quantities (0p, /3), My {(j}p, /3) and 

Mry {'Pp: /?) = 4f \J ~ ^^^^ ^® functions of (/ip and /3. Contour 

plots of these three functions are shown in Figure |5] Again, agreements and 
disagreements in symmetries are found by comparison with Figure [T] Inspec- 
tion of Figures|4]and|5]immediately reveals differences in the transformation 
behaviour of p^ and p°^. 
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As shown earlier it is possible to implement any logic gate from a NMR 
measurable quantity that has the symmetry of a class 3 gate in Figure [T](e.g. 
Mx ((f>p,P) in Figure |4ja)). However, none of the three contour plots in 
Figure [5] displays this symmetry and no class 3 gates can be implemented. 
The symmetries present belong to those of class 0, 1 and 2 gates and can be 
found in any of the three contour plots. Implementing class 1 gates requires 
taking advantage of the periodicity of the functions {4>p, /3), My (0p, (3) 
and Mxy {(pp, f3) in order to generate canalising inputs. Class 2 gates use a 
reduced parameter range for (j)p and /3. 

Permutation of (f>p and /3 changes {ipp, f3), My (0p, (3), while M^y {4>p, /3) 
stays invariant under this operation. M^ {(pp, /3) and M^y {(j)p, (3) are invari- 
ant under input inversion 0p — i/ip or /? — > — /3 while My {(j)p, (3) is only 
invariant under /3 ~^ —/3. All functions are now periodic with tt in 0p be- 
cause of their dependence on sin^ (j)p or sin 2(/)p. AI^. {(j)p, (3) and My (0p, (3) 
have a 2it periodicity interval in /3 and tt periodicity for M^y (cfip, (3). All 
the functions M^ {(pp, f3). My {(pp, (3) and M^y {(pp, (3) are the products of 
two trigonometric functions of which at least one is quadratic. The squaring 
reduces the trigonometric function values to the interval [0, 1] and therefore, it 
is not possible to have a symmetry behaviour required for class 3 gates. 

The rotation operator iJ^^ (/?) not commuting with the initial density matrix 
p'^^ leads to the functions M^ {pp, (3), My {(pp, j3) and M^y {(pp, (3) described 
by products of linear trigonometric functions and the ability to implement 
all logic gates. The rotation operator i?^^ {f3) commuting with the initial den- 
sity matrix p^, however, leads to the functions M^ {(pp, (3), My {(pp, (3) and 
Mxy {(pp, (3) characterised by quadratic trigonometric factors and does not per- 
mit the implementation of all logic gates. Achieving universality when starting 
from px therefore requires a different rotation operator that never commutes 
with px ■ This could be achieved by choosing different pairs of parameters 
from the set of parameters describing the single-pulse experiment (Figure |3]l 
such that a non-commuting rotation operator for p^ results. Alternatively, one 
can keep the pair of parameters (pp and j3, and use more than one r.f. pulse, 
with rotation operators such as i?0p2 {(32) ^^pi {Pi)- This latter option we 
consider next. 

Two-pulse gates 

The initial density matrix p^. transforms under the influence of two r.f. pulses 

i?0p. (/32)i?0p, (^i)as 

P^^^ (0p2, /32, 0pl, A) = {P2) i?0p, P.Rl^, {Pi) {132) (29) 
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Thus, the NMR detectable quantities {4>p2, f32, <t>pi, f3i). 
My (</'p2,/32,(/)pi,/3i) and M^y {(f>p2,(32,4'pi,l3i) are functions of the four 
variables 4>p2, P2, ^pi and /3i. Any two of these variables can serve as logic 
gate input parameters with the remaining two fixed. This offers a high degree 
of flexibility in assigning parameters as variable (and to control the logic gate), 
and as fixed. 

A trivial example is based on the fact that all 2-input gates can be imple- 
mented when starting from thermal equilibrium magnetisation by a single 
pulse. Starting from p^., a first r.f. pulse Ril (— f ) generates pf^ and yields 
0pi = ^ and /3i = — ^. From there on any second r.f. pulse R(jj^^ {P2) gen- 
erates a functional behaviour identical to that of the single-pulse scenario in 
Eq. ( |25] l, generating universality. 

Another strategy uses the commutation properties of the rotation opera- 
tors. In the single-pulse scenario, starting from p°^ all logic gates can be 
implemented because the overall rotation operator never commutes with the 
initial density matrix. In the two-pulse scenario, i?<^p2 {h) -R^pi {Pi) must 
not commute with p^- This implies further that the first r.f. pulse Rcj,-^^ (/?i) 
must not commute with px, otherwise Eq. (|29| simplifies to p*^^' {4>p2Th) = 
R<l}p2 PxR^p 2 {P2) (identical to Eq. ( p7j l) which can not implement all 
logic gates. To avoid commutation, (/ipi should never be fixed to values mr 
(Eq. (|26|). The second r.f. pulse Rij,^^ (/^s) must not commute with the spin 
basis operators or ly when calculating {(f>p2, P2, 0pi, Pi) = (Ix) or 
My ((/)p2, /32, 0pi, /3i) = {ly) (Eq. (|23t), otherwise 

(/32) -R^pi (/3i) Px-R^pi Wi) -fl^p2 W2) is just a similarity transform 
of i?0pi iPi) PxRIp J iPi) Ix under which the trace is invariant. A com- 
mutation [R^p2 iP'i) '^2;] =0 can be avoided by setting (f)p2 7^ nir. For 
My ((/)p2, /32, 0pi, /3i) a commutator [R^,^^ {h) ,Iy\ =0 is avoided for (/)p ^ 
(2n + 1) I . These constraints on the parameters of type phase 0p imple- 
ment the functions {4>p2^ /32, '/'pi, such as those shown in Table|5]and 
Figure |6] these examples were deliberately chosen such that the two fixed 
variables are taken as having equal values of |. 

One can see that class 0, 1 and 3 gates can be implemented by all possible 
permutations of fixed pairs of variables of identical value (here The 
non-commutation constraint provides the non-canalising behaviour of class 3 
gates, while the inherent periodicity of the system is sufficient to provide the 
canalising input configurations necessary for class and 1 gates. However, 
periodicity does not suffice for a class 2 gate implementation. Comparing 
the single-pulse scenario (Figure |4]^a)) with the contour plots in Figure |6] 
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7r 
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— ^ sint/ipi sin/3i 



TABLE 5 

Examples of logic gate implementations from Mx {4'p2 , /32 , 0pi , /3i ) using any possi- 
ble pair of variables as gate input A and B, and fixing the remaining two parameters 
to ^ . Apart from a change in sign, examples (e) and (f) are identical to the case in 
Figure|4|a). Example (d) is identical to (b). 



(b) (c) 




0p2[rad]/Tr 0p2 [rad]/Ti ^3 [rad]/7T 



FIGURE 6 

Contour plots of {4>p2, P2, 4>pi, /3i) corresponding to examples (a)-(c) in Tablejs] 
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(a) (b) 




4ip2[rad]/n ifp2[rad]/7i 



FIGURE 7 

Contour plot of {<i>p2,l32, <i>pi,l3i) for (a) /3i = f , /32 = vr and for (b) (?!>pi = f , 
02 =n. 



demonstrates that class 2 gate implementations require horizontal and vertical 
traces of constant value zero. Such traces are absent in the contour plots in 
Figure |6j In short, arbitrarily choosing pairs of fixed parameters as equal 
valued (not just for ^) leads to the loss of universality. 

Let us examine if universality is regained if the fixed parameters are no 
longer taken as equal valued. For example, taking /3i = ^ and /32 = tt gives 
{(j)p2,P2, 0pi, = ^ COS (20p2 - (j)pi) COS (0pi). The corresponding 
contour plot is shown in Figure |7]^a) highlighting the presence of constant 
horizontal zero-valued traces. Allowing for unequal-valued pairs of parameters 
/3i and /32 (or (ppi and ) permits implementation of class 0, 1 and 2 gates 
but not class 3. 

Universality is regained if pairs of parameters of different types, e.g. (ppi and 
f32 are fixed and do not assume equal values. Figure|7j;b) shows that 0pi = ^, 
f32 = TT gives Afj. (0p2,/32,0pi,/3i) = ^ cos 20P2 cos /3i . Both horizontal 
and vertical zero-valued traces are re-established in the corresponding contour 
plot (Figure [7jb)); all 2-input logic gates can be obtained. 

The flexibility and the ease of implementation as seen from our illustrative 
examples are good indicators for the richness of the natural computational 
potential of this system. Computational operations more complicated than just 
2-input boolean logics are therefore well within the grasp of this system. This 
includes multi-input logic gates, continuous logic, and analogue computing 
implementations. The ubiquitous occurrence of function values {—1, 0, 1}, 
due to the trigonometric system functions, especially holds promise for the 
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implementation of a balanced ternary logic as being natural to this particular 
system |7]. Note also, that we can reverse the operation of a logic gate if we 
know the pulse(s) that originally generated the output [20 1. This does not 
match with the usual predictions about a computation where e.g. NAND gates 
are not reversible |8^|. Obviously our raw output stores additional information 
that can be exploited for more sophisticated computations. 

4 SUMMARY AND CONCLUSIONS 

We have provided a design approach to analysing novel substrates in order 
to determine which of their parameters can be used to implement boolean 
logic gates. We have illustrated this with a case study drawn from NMR-based 
classical computation. The design process requires cataloguing the parameters 
that are naturally used to describe and manipulate the target system, analysing 
their behaviours in combinations, and then matching the resulting patterns of 
behaviour with the corresponding behaviour patterns of the target gate classes. 
Our case study here focusses on two particular parameters {(3 = KpTp and 
(pp); a full design study would assess other combinations of other parameters, 
since the aim is not simply to find some solution, but to analyse the 'natural' 
computational capabilities of the substrate. For example, in our work on 
continuous gates [6J, we focussed on Wp and Tp. 

The design approach is not restricted to 2-input gates; combining several 
parameters can produce patterns corresponding to more complex gates. More 
sophisticated experiments could directly exploit symmetry properties of the 
Hamiltonian (Eq. ([TT|-([T4l)). Eventually, however, specific behaviours will be 
more easily achieved by combining simpler gates in circuits than by directly 
designing complex gates. 

Circuit design requires additional analysis to determine how individual 
gates can be combined in a circuit in a manner 'natural' for the substrate 
in question. Circuit design requires determining techniques for: sequencing 
parameter manipulations to implement the sequence of gate operations in a 
circuit; combining manipulations to implement multiple gate operations in 
parallel; routing and transforming the output of one gate to the appropriate 
input of the next. We will address these circuit design issues in a companion 
paper f5]. 

This work addresses only how to implement 'classical' boolean logic gates 
in unconventional substrates, yet it is clear that the example system has further 
information available, that is thrown away when viewing it as a boolean gate. 
If this extra information is instead retained and exploited, more powerful 
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computation becomes available fV\. The real computational power of these 
novel substrates will come from not viewing them as merely altemative ways 
of implementing classical logic gates, but from exploiting them to implement 
non-classical forms of computation: quantum, continuous (analogue), hybrid, 
and more. The design approach described here forms a first step in a principled 
approach for analysing substrates with a view to performing a specified form 
of computation. 
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